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ABSTRACT 



Aims. A Gaussianity analysis using a goodness-of-fit test and the Minkowski functionals on the sphere has been per- 
formed to study the measured Archeops Cosmic Microwave Background (CMB) temperature anisotropy data for a 143 
GHz Archeops bolometer. We consider large angular scales, greater than 1.8 degrees, and a large fraction of the North 
Galactic hemisphere, around 16%, with a galactic latitude b > 15 degrees. 

Methods. The considered goodness-of-fit test, first proposed by Rayner fc Best (1989), has been applied to the data 
after a signal-to-noise decomposition. The three Minkowski functionals on the sphere have been used to construct a % 2 
statistic using different thresholds. The first method has been calibrated using simulations of Archeops data containing 
the CMB signal and instrumental noise in order to check its asymptotic convergence. Two kind of maps produced with 
two different map- making techniques (coaddition and Mirage) have been analysed. 

Results. Archeops maps for both Mirage and coaddition map-making, have been found to be compatible with 
Gaussianity. From these results we can exclude a dust and atmospheric contamination larger than 7.8% (90% CL). 
Also the non-linear coupling parameter f n i can be constrained to be f n i = 200 lgJo° a ^ the 95% CL and on angular 
scales of 1.8 degrees. For comparison, the same method has been applied to data from the NASA WMAP satellite in 
the same region of sky. The 1-year and 3-year releases have been used. Results are compatible with those obtained with 
Archeops, implying in particular an upper limit for f n i on degree angular scales. 

Key words. Cosmology - data analysis - observations - cosmic microwave background 



1. Introduction 

k/ ' Accordin g to the inflationar y un i verse theory (s ee for 
' example iGuthl Il98ll : iLindd . Il990l : iLyth fc Riottol . Il998l : 
■ iLiddle fc LvthL 120001 ), the primordial density fluctuations 
$-h ' are distributed following very precisely a Gaussian proba- 
5^ , bility density function (pdf). These fluctuations in the mat- 
ter density will produce anisotropics in the temperature of 
the Cosmic Microwave Background (CMB) whose pdf is 
also Gaussian. In this manner, when the Gaussianity of the 
CMB radiation is analysed the standard inflationary the- 
ory is tested as well as its alternatives (for example cosmic 
strings) which generically predict deviations from it in dif- 
ferent ways. In addition, the search for non-Gaussianities 
has become a powerful tool to detect the presence of resid- 
ual foregrounds, secondary anisotropics (such as gravita- 
tional lensing, Sunyaev-Zel'dovich effect) and unidentified 
systematic errors, which leave clearly non- Gaussian im- 
prints on the CMB-anisotropies data. There are many tech- 
niques to test Gaussianity, many of them developed previ- 
ously as general statistical methods to test the normality of 
a data set, and others specifically for the CMB anisotropics. 

Among those methods, we can mention the esti- 
mator for non- Gaussian it y ba se d on th e CMB bis- 
pectrum ([Ferreira et all Il998l : iMagueiiol, 120001), ge- 
ometrical estimators on the sphere (Barreiro et all 



120011 iMonteserfn et all I2QQ5L l2QQ6h M inkowski function- 
als ([Gott et all ll99Ct IKomatsu et all 12003th gq o dness- 



of-fit tests (IRavner fc BestL Il989[ ~ Aliaga et al.1 , 
Barreiro et all |2006), wavelets (iFerreira et al 



2003; 



1997 



Hob son et all ll999l ; lBarreiro et all l2QQQh and steerable fil- 



ters to search alignment structures (jWiaux et all 120051 ) . 

Some of them have been applied to the CMB providing 
different results. For example WMAP data are compati- 
ble with Gauss i anity according to the W MAP team (see 
IKomatsu et~afl 120031 : ISpergel et all . l2Q07h whereas others 
have found eviden ces of non-Gauss s ianiti es in the same 
WMAP maps, like ICopi et all (|2004 l2006f ) (using a tech- 
nique called multipole vector framework), lEriksen et al.l 
(|2004 l2005f ) (finding asy mmetries using local es timators of 
the n- point correlations), IVielva et al.l (|2004t ) ; ICruz et all 
(I2QQ5L I2QQ6L |2QQ7h ( the Cold Spot detected with wavelets), 
lLarson et al.l (120041 ) (cold and hot spots different from 
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the ones expected in Gaussian temperature fluctuations), 
among others. 

In thi s work the smooth go odness-of-fit test first pro- 
posed by IRavner fc Best! (|l989f ) (hereafter R&BT) will be 
implemented to analyse the Gaussianity of the Archeops 
data. This method has been al r eady a pplied successfully 
to the MAXIMA (ICayon et all l2Q03bh and VSA exper- 
iments ([Aliaga et all I2QQ5I : iRubino-Martfn et all l2006f ). 
The Archeops data will be as well analysed with the 
morphological descriptors known as Minkowski functionals 
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([Schmalzing fc Gorskl [l998l : iGott etallll990h - The idea is 
to use both methods in the Gaussianity analysis for com- 
parison of the sensitivities of the two techniques and cross- 
checking of the results on the amount of dust contamination 
and the amplitude of the non-linear coupling parameter. 

This is the first analysis of Gaussianity of the Archeops 
experiment data. We have analysed the data for one of 
the Archeops bolometer at 143 GHz. This bolometer is the 
most sensitive and one of the most relevant for CMB ob- 
servations. As a complementary analysis, we present the 
results of the same goodness-of-fit test applied to WMAP 
data with approximately the same mask as the one used for 
Archeops. The purpose is to check whether the results are 
consistent for both data sets. 

This paper has the following layout: in Section 2 
the R&BT applied to signal-to-noise eigenmodes and the 
Minkowski functionals are described. The experiment, main 
properties of data sets and masks are summarized in Section 
3. Section 4 is dedicated to the calibration and checking 
of both methods with some "realistic" CMB anisotropy 
Gaussian simulations, where we know in advance the out- 
put of the techniques. Section 5 contains the Archeops data 
analysis as well as results. In Section 6 WMAP 1-year and 
3-year data are analysed and compared with Archeops re- 
sults. Finally in Section 7 the main conclusions are pre- 
sented. 



2. Goodness-of-fit tests and Minkowski functionals 

In this section, on the one hand, we will describe briefly the 
"goodness-of-fit technique" applied to test the Gaussianity 
of a set of "signal-to-noise eigenmodes" derived from mea- 
surements of the CMB temperature anisotropics. On the 
other hand, we will explain the Gaussianity analysis based 
on the Minkowski functionals. 



2.1. Smooth tests of goodness-of-fit 

Given a set of n random numbers, {yi}\=i, it is sometimes 
interesting to check whether they behave statistically ac- 
cording to one specific pdf, f(y,0), that is, if the probabil- 
ity of finding a random number y in an interval between 
y and y + Ay, with Ay > 0, is given by f(y o ,0)Ay. A 
scalar or vector variable is introduced, which allows us to 
move smoothly between different pdf 's in their correspond- 
ing space of normalized functions. 

This statistical analysis consists in testing the null hy- 
pothesis, Ho : {0 = 0} against the alternative hypothesis, 
/\' : {« — ()}. 

Inside the family of smooth goodness-of-fit tests, we can 
consider an order k alter native pdf Qkiv charact erized 
by a pdf of the form fsee lRavner fc Best! . I1989L ll99Qh 



g k (y,0)=C(0)exp 



i=i 



f(y) 



(i) 



is a set of k parameters to smoothly cover our space of 
pdf's, f(y) is the null hypothesis pdf (e. g. the Gaussian dis- 
tribution), hi(y) form a complete set of orthonormal func- 
tiomQ on f(y) and C{6) is a normalization constant. 



The "score statistic" is used to evaluate the simple null 
hypothesis Hq. With this statistic one can estimate the sta- 
tistical significance of through the " Maximum Likelihoo d 
Method". Following the notation bv lAliaga et all (|2003f h 
the score statistic for this goodness-of-fit test is 



and the Uf quantities are given by 



h i(Vj) 



(2) 



(3) 



In the case of a Gaussian pdf, hi(x) are the "normalized 
Hermite-Chebyshev polynomials" . If the null hypothesis is 
satisfied then the Ui quantities have a statistically normal 
behaviour and therefore Uf behave like a Xi distribution 



(4) 



It is possible to write the Uf statistical quantities in terms 
of the moments of order k derived from the set of n ran- 
dom numbers to be analysed, = 1/^ J^? =1 2/? , (see for 



example lAliaga et all I2QQ3L 120051 ) . 

In this work, the five first statistics Uf have been used 
and can be related to the /c-order moments in the following 

way, 



ul 
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(5) 



1 I-oo h i(y)f(y) h j(y) d y = h 



The first few statistics are generally the most sensitive for 
most of the applications. In our case higher order Uf statis- 
tics are dominated by errors (because of usual propagation 
of errors) and therefore are not very useful in practice. This 
will be described in detail in section [H 

2.2. Signal-to-noise eigenmode analysis 

At this point, we have described the method that will be 
used to analyse a set of n random numbers to test whether 
their pdf is the normal distribution or not. 

The next step is to compute the set of numbers to be 
analysed. In our case they come from the so called "signal- 
to- noise eigenm odes" , firstly introduced in the CMB field 
bv lBondl (pL995). Our observational data, (the fluctuation in 
the temperature of the incoming blackbody radiation mea- 
sured for each direction n in the sky, AT(n)/T), can be 
interpreted as originated from several sources: all the emis- 
sions coming from the sky (CMB signal, Galactic and ex- 
tragalactic foregrounds and a tmosphere) and the measu red 
instrumental Gaussian noise (Maci as- Perez et all 12007ft . 

The total area observed by the experiment is usually 
divided in equal area pixels identified by their center di- 
rection n and to which the measurements, AT(n)/T, are 
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assigned. To obtain the "signal-to-noise eigenmodes" , we 
expand the pixel values of the map, AT(n)/T, into a linear 
combination in which the transformed instrumental noise 
(hereafter the noise) and the transformed theoretical CMB 
signal (hereafter the signal) do not have correlations. 

For the signal-to-noise decomposition it is necessary to 
calculate signal and noise covariance matrices. The temper- 
ature covariance between two pixels i and j is given by 

dj = {ATiATj) - {AT^iATj) (6) 

where the brackets () represent the average over several re- 
alizations of temperature anisotropy maps. Thus we can 
construct the signal (noise) covariance matrices, S (iV), 
averaging on signal AT s (n) (noise AT n (n)) realizations. 
Since the data represent temperature fluctuations around 
the mean then it is trivially satisfied that (AT s (n)) = 
(AT n (n)> = 0. Therefore, dj = (ATiATj), the correlation 
matrix. 

Once we select a set of n directions in the sky (pixels) 
and we construct S and N matrices, which have the same 
dimension n and are symmetrical, we can compute the so 
called "signal-to-noise matrix" A 

A = L-JSL-J (7) 

where is the Cholesky matrix of N, defined as N = 
LnL 1 n . Ln can be obtained from the diagonalization of the 
N matrix. Suppose Dn is the diagonal matrix of eigenval- 
ues of N, and Rn a matrix of the eigenvectors of N, related 
by R n NRn = Dn. Then it is satisfied that Ln = RnD]^ 2 , 

1 /2 

where is the square root matrix of Dn- 

If d is the vector of dimension n representing the data 
assigned to the pixels in the sky, the signal-to-noise eigen- 
modes can be written as 

£ = E^Ljd (8) 



2.3. Minkowski functional 

Considering the temperature anisotropics of the CMB as 
a scalar field on the sphere we can define the set of 
coordinates Q v where AT(n) > v for a given thresh- 
old v, and its complement ary set V v . As it is stated in 
ISchmalzing fc Gorskil (|l998l ). any morphological descriptor 
on the sphere is a linear combination of 3 Minkowski func- 
tional. These functionals are: the area A(y) of the excur- 
sion set Q V1 the contour length C{y) of the excursion set 
Qzy, and the genus G(y) (defined as the number of hot spots 
above v minus the number of cold spots below that thresh- 
old). 

For a Gaussian random field the mean values of these 
functionals are 



(A{y)) = 

(C(u)) = 
(G(u)) = 




(10) 

where r is a parame t er related with the coheren ce angle 
(jBarreiro et all 1200 it ISchmalzing fc Gorskil . Il998h . 

The Gaussianity test with the Minkowski functionals 
is performed through a y 2 test as described for example 
in lKomatsu et all (|2003t ): ISpergel et al.l(|2007l ). Considering 
n t h possible thresholds v we can define a 3n t h vector v = 
(A{y),C(y),G(y)). The x 2 statistic is then defined 

X 2 = 5>(i) - (v(i)))C^(vV) - (v(j))) (11) 

where (v(j)) is the expected value of v(j) and C is the 
corresponding covariance matrix for all possible thresholds 
and functionals. 



where Ra is the matrix of eigenvectors of A and Da the 
diagonal matrix of eigenvalues of A, R\ARa = Da- 

The yi quantities to be analysed with the goodness-of-fit 
test defined in the previous section are 



Vi 



V 1 + {PA)i 



(9) 



It can be easily demonstrated that if the vector of data d 
satisfies (d) = then (yi) = 0. In the case AT = AT S + 
AT n , from the definition of signal-to- noise eigenmodes in 
(|8j), the definition of yi in ([9]), and properties of correlation 
matrices it follows that (yf) = 1. 

Supposing that the original map d is multi-normal, then 
our {yi} numbers keep the Gaussian character because both 
set of numbers are connected by linear operations. More 
precisely, they follow a normal pdf with zero mean and unit 
variance, iV(0, 1). Moreover, for different indexes i and j, 
yi and yj are independent. 

Finally, for Gaussian data d each Uf statistics, defined 
in (|3J), is distributed as a xl- The decision to accept or reject 
the null hypothesis will be therefore based on this pdf, as 
will be seen in sections 5 and 6 when the test is applied to 
the Archeops and WMAP data. 



3. Archeops data sets 

3.1. The Archeops experiment 

Archeop^3 is a balloon borne experiment dedicated to mea- 
sure the CMB temperature anisotropics from large to sma ll 
angular scales (Ben oit et all l2QQ3at iTristram et ail . [2005). 
It has given the first link in the Ci det ermination between 
the COBE large angular scales data (|Smoot et all Il992f ) 
to the first acou stic peak as measu r ed by BOOMERanG 
and MAXIMA (|de Bernardis et all 120001 : lHananv et "all 
l2Q00h . Archeops was also designed as a test bed for the 
forthcoming Planck H igh Frequency Instrument (HFI), 
([Lamar re et all l2QQ3h . Therefore, Archeops shared with 
Planck the same technological design: a Gregorian off-axis 
telescope with a 1.5 m primary mirror, bolometers oper- 
ating at 143, 217, 353 and 545 GHz cooled down at 100 
mK by a 3 He/ 4 He dilution designed to work at zero grav- 
ity and a similar scanning strategy. Archeops was launched 
on February 7 th , 2002, from the CNES/Swedish facility 
of Esrange, near Kir una (Sweden). 12 hours of high qual- 
ity night data were gathered. This data corresponds to a 
coverage of approximately 30% of the sky, including the 
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Galactic plane. More details about t he instrument and th e 
flight performance can be found in iBenoit etafl ([2003b); 
IMacfas-Perez et all (|2007l ). From its four frequency bands 
the two lowest (143 and 217 GHz) were dedicated to the 
observation of the CMB and the others (353 and 545 GHz) 
to the monitoring and calibration of both atmospheric and 
Galactic emissions. 

In the following, we focus on the analysis of the most 
sensitive 143 GHz Archeops bolometer that also presents 
the lowest level of contamination by systematic effects. 

Although the Archeops resolution is typically of 
10 arcmin, for this analysis we are interested in the 
Gaussianity of the large angular scale anisotropics. 
Therefore, we decided to u se low resolution maps at 
HEALPix ( Go rski et al.Ll2QQ5h N s ^ e = 32 to consider scales 
above 1.8 degrees. 

3.2. Data processing 

We describe here briefly the way that Archeops data 
were processed. For a more detailed description see 
IMacfas-Perez et al.l ([20071 ). 

In the Time Ordered Information (TOI) corrupted data 
are flagged (representing less than 1.5% of the whole data 
set). Low frequency drifts, correlated to house-keeping data 
are removed using the latter as templates. A high frequency 
decorrelation is also performed to remove some bursts of 
non- stationary high-frequency noise. Corrected timelines 
are then deconvolved from the bolometer time constant and 
the flagged corrupted data are replaced by a realization of 
noise. Finally, low time frequency atmospheric residuals are 
subtracted using a destriping procedure which slightly fil- 
ters out the sky signal to a maximum of 5%. 

Archeops cleaned TOIs at 143 GHz are contaminated by 
atmospheric and Galactic dust residuals, even at interme- 
diate Galactic latitudes. Atmospheric residuals contributes 
mainly at frequencies lower than 2 Hz in the timeline and 
follows approximatively a v 2 law in antenna temperature. 
Galactic dust presents a grey body spectrum at about 17 K 
and with an emissivity of about v 2 . To suppress both resid- 
ual dust and atmospheric signals, data are decorrelated us- 
ing a linear combination of the high frequency photometric 
pixels (353 and 545 GHz) and of synthetic dust timelines. 

We have used in this work two kind of map-making for 
the TOIs of Archeops data and of the simulations. The 
first one is an optimal map- making procedure called Mirage 
([Yvon fc Mavetll2005h . Mirage is based on a two-phase it- 
erative algorithm, involving optimal map-making together 
with low frequency drift removal and Butterworth high-pass 
filtering. A conjugate gradient method is used for resolving 
the linear system. The second is a procedure that performs 
coaddition. This means that all the TOI points correspond- 
ing to a given pixel are averaged. 

To produce a CMB simulation, a random CMB map 
with the power spe ctrum of the Archeops model (see 
Ben oit et al.l ([200 3b) and figure [T]) is generated and from 
this map an Archeops TOI is produced. This TOI is treated 
with the two map-making methods described above to pro- 
duce a map. To perform a noise simulation we produce 
a Gaussian constrained realization of the Archeops noise 
power spectrum in the time domain. The TOI produced 
this way is then projected into a map using the above map- 
making techniques. 



-i • I | B< t Fit P< i | tru 




Fig. 1. Archeops Best Fit Power Spectrum used to simulate 
the Archeops CMB signal. 
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Fig. 2. Mirage Archeops data from the best bolometer at 
143 GHz presented at HEALPix resolution N 3 id e = 32, (« 
1.8 degrees). This map is centered on Galactic longitude 
I = 180 degrees. Galactic and South Equator pixels have 
been masked. Grid lines are spaced by 20 degrees. 

The analysis has been performed on a fraction of the 
Archeops observed region masking out pixels with Galactic 
latitude below 15 degrees, \b\ < 15°. The southern sky data 
were not included in the analysis as they are more contam- 
inated by systematics in the form of residual stripes com- 
ing from the Fourier filtering and destriping o f the data 
in the time domain ([Macfas-Perez et all 120071 ) that pro- 
duces ringing around the Galactic plane. In the case of the 
CMB power spectrum analysis presented in lTristram et al.l 
([20051 ) this southern sky region was used as increased sig- 
nificantly the signal to noise ratio at small angular scales 
which are not affected by this systematic effect. This is not 
the case for the analysis presented in this paper where we 
are more interested in large angular scales where this sys- 
tematic becomes important. In figure [2] we plot the region 
of data considered for the analysis. These data correspond 
to 1995 pixels (16% of the sky) from a total of 12288 pixels 
for a complete map at this resolution. 

4. Calibrating the method: analysis on Gaussian 
simulations 

To develop the R&BT non- Gaussianity test, it is necessary 
to calculate the signal (S) and the noise (N) correlation ma- 
trices among the selected pixels. We computed these ma- 
trices averaging on simulations by means of equation [6l For 
this purpose Monte Carlo Gaussian simulations of Archeops 
CMB signal and instrumental noise were produced. The 
number of performed simulations for the map generated 
with the Mirage map-making procedure were 2.86 x 10 5 
for the signal and 2.75 x 10 5 for the noise, whereas for the 
coaddition procedure they were 5 x 10 5 and 5 x 10 5 for the 
signal and noise respectively. 90 dual-core 3.2 GHz proces- 
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Fig. 4. From left to right, mean and dispersion of Uf statistics (where i goes from 1 to 5) for different signal-to-noise 
cuts, corresponding to 10 4 signal plus noise Mirage simulations. 




Fig. 5. From left to right, mean and dispersion of Uf statistics (where i goes from 1 to 5) for different signal-to-noise 
cuts, corresponding to 10 4 signal plus noise coaddition simulations. 




Fig. 6. From left to right, and from top to bottom, distribution of the Uf statistics, from a set of 10 4 Gaussian 
Mirage simulations analysed in the same region than the data (figure [2]). The signal-to- noise cut which has been used is 
(s/n) c = 0.30. Solid lines are the theoretical distribution (xi) normalized to the number of simulations and the size of 
the binned cell. 



sors from the IFCA computing facilities were used. Each 
Mirage simulation took 180 s of real CPU time and 1.0 GB 
of RAM memory, whereas these values were 70 s and 0.04 
GB respectively for each coaddition simulation. 

The high number of simulations and the corresponding 
computational requirements were needed to achieve conver- 
gence in the construction of the correlation matrices. The 
main reason for the low convergence relies on the specific 
properties of our correlation matrices. Archeops noise is 



correlated at large scales, which means that the TV matrix 
is neither diagonal nor sparse. The Archeops signal corre- 
lation matrix contains correlations at large scales for which 
the convergence is much slower than for the small scales 
due to the cosmic variance. In both cases many simulations 
(~ 10 5 ) were required in order to compute these matrices. 

One way to quantify the degree of convergence of these 
matrices is by analysing Gaussian simulations. The Uf 
statistics for a set of Gaussian simulations should have a 
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Fig. 7. From left to right, and from top to bottom, distribution of the Uf statistics, from a set of 10 4 Gaussian 
coaddition simulations analysed in the same region than the data (figure [2]). The signal-to- noise cut which has been used 
is (s/n) c = 0.27. Solid lines are the theoretical distribution (xD normalized to the number of simulations and the size of 
the binned cell. 



Area Contour Length Genus 




Fig. 8. From left to right, mean values of the three Minkowski functionals and their corresponding error-bars for a 
set of 1000 (noiseless) CMB Gaussian simulations. These simulations have been generated using Archeops best fit power 
spectrum and have not been masked. Notice the good agreement between the theoretical predictions and the results 
obtained from simulations. 




Fig. 3. Number of normalized signal-to- noise eigenmodes 
Ui for which their associated A matrix eigenvalues, (Da)i, 
satisfy (D A )i > (s/n) 2 c . 



Xi pdf. This can be tested, for example, by calculating the 
mean and the variance of the Uf statistics for 10 4 Gaussian 
signal plus noise simulations. For the Gaussian case, the 



mean should be equal to 1 and the dispersion equal to y/2 
(this is the null hypothesis, Hq). 

Following lAliaga et all (|2005f ): iRubiho-Martm et all 

(2006), the Uf are computed for a subset of signal-to- noise 
eigenmodes which are those associated with eigenvalues of 
the signal-to-noise matrix A satisfying (Da)% > (s/n)^ 
where (s/n) c is a given signal-to-noise ratio cut. In figure 
[3] the number of eigenmodes {^}'s, which obey {Da)% > 
(s/n)l, in terms of s/n is plotted. 

In figure 0] we show the mean and dispersion of the five 
first Uf statistics for different signal-to-noise cuts corre- 
sponding to all possible eigenvalues of the A matrix. The 
Uf values come from a set of 10 4 Gaussian Archeops sig- 
nal plus noise Mirage simulations. It can be seen that for 
( s / n )c ^ 2 mean values are close to 1 and the dispersion 
close to \]~2 (except for the f/f statisti c whose dispersion i s 
always larger than 2). As shown by e.g. lAliaga et al.l (|2QQ5h . 
the expected value of Uf is equal to 1 independently of the 
number of {yi} used. This explains why we have got the 
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mean of Uf very close to 1 for every signal-to-noise cut. The 
dispersion is equal to y/2 asymptotically, when the number 
of {yi} used is high. In our case, this happens for low signal- 
to-noise cuts, when enough {^}'s are used to compute the 
statistics. In figure [5] the same quantities have been plotted 
for the 10 4 Gaussian Archeops signal plus noise coaddi- 
tion simulations. Similar conclusions can be derived in this 
case. Notice, however, that the results are closer to theoret- 
ical values when the analysis is performed using the Mirage 
maps. In this case the correlation matrices have converged 
with less simulations than in the coaddition case. This is 
one of the advantages of using Mirage simulations over the 
coaddition ones, although the production of a Mirage map 
requires more CPU time and RAM memory than a coaddi- 
tion map. 

Since the computation of high order Uf statistics in- 
volve high powers of the eigenmodes, the convergence of 
their dispersion to the theoretical values at a given (s/n) c 
is slower than for the low order ones (as can be seen in the 
right panels of figures [4] and [5]) . 

A more exhaustive check for the convergence of the Uf 
statistics is done by comparing their theoretical pdf with 
the histograms obtained from the simulated data. Given a 
signal-to-noise ratio cut (s/n) c for the calculation of the Uf 
statistics, it is possible to make a histogram with the cor- 
responding values of the Uf statistics from the same sets 
of 10 4 simulations. Figure |B] compares the histograms for 
the first five statistics calculated using all the eigenmodes 
(s/n > 0.30) for the Mirage simulations with the theoretical 
expectation of a xi distribution. In table [1] the mean and 
the dispersion of these histograms are presented. In figure 
the same comparison is shown for the coaddition simula- 
tions considering also all the eigenmodes (s/n > 0.27). The 
corresponding mean and dispersion of these histograms are 
given in table El 

In summary, the four statistics U%, U$ , and U% have 
a pdf's compatible with the theoretical one whereas C/f 
starts to deviate from it. The discrepancy, already present 
in the dispersion, increases for higher orders. The reason is 
that high order moments enlarge possible errors present in 
the computed correlation matrices and are propagated in 
the diagonalization processes. In any case, the f7f statistic 
can still be used for the Gaussian analysis if the distribu- 
tion obtained from the simulations, instead of the theoret- 
ical one, is used as reference. Although this is not as opti- 
mal as using the theoretical xi distribution, it is however 
a good compromise taking into account the huge computa- 
tional resources needed to produce a very large number of 
simulations. 

For the Minkowski functionals analysis the expected val- 
ues given by equation [10] cannot be applied to our problem 
because of the contour restrictions of the mask and the 
presence of anisotropic noise. Nevertheless in order to test 
our Minkowski functional codes we performed an analysis 
on (noiseless) CMB Gaussian simulations over all sky and 
1.8 degrees resolution generated using the best fit Archeops 
power spectrum. Analysing them for thresholds from —2.5a 
to 2.5a (where a is the standard deviation of the corre- 
sponding simulation), we obtained that the results from 
simulations are compatible with the theoretical predictions 
(see fig. E]). 



Table 1. Mean and dispersion of Uf statistics from 10 4 
Mirage simulations for a signal-to- noise ratio cut of 0.30. 
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Table 2. Mean and dispersion of Uf statistics from 10 4 
coaddition simulations for a signal-to- noise ratio cut of 0.27. 
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5. Gaussianity test on Archeops data 

We have applied the R&BT to the Archeops 143if03 
bolometer map. The signal-to-noise eigenmodes have been 
computed with the correlation matrices described in section 
IH for each map- making case. We have checked in that sec- 
tion that these signal and noise matrices provide Uf statis- 
tics compatible with Gaussianity for Gaussian simulations. 

We have applied this test to the Archeops data for 
the Mirage and coaddition map-making. The Uf statistics, 
computed for the 1995 pixels of the previously described 
Archeops data, are displayed in figures [9] and [10J The Uf 
statistics are plotted, from i = 1 to 5, versus the signal-to- 
noise eigenmode cut. 

For the Mirage map-making, results are displayed on 
figure [9l We can see that all the Uf statistics are below 5 
for all the signal-to-noise cuts. This means that the data is 
compatible with Gaussianity. 

For coaddition map-making, we can see from figure 
[10] that whatever the signal-to-noise eigenmode cut is, Uf 
statistics for the 143K03 bolometer data are below 5, ex- 
cept for Ul for signal-to- noise cuts below 0.5. It reaches the 
maximum value of 7.97 at the minimum signal-to-noise cut 
of 0.27. The upper tail probability^ for C/| = 7.97 from the 
Xi distribution (equation [4]) is 0.5%. Comparing with the 
set of coaddition Gausian simulations we found that this 
upper tail probability is 0.6%, (see table [3]), in good agree- 
ment with the theoretical expectation. Nevertheless, as we 
have computed Uf statistics for all possible signal-to-noise 
cuts, it is important to estimate the significance of finding 
any simulation with Ul > 7.97 in at least one of them. This 
is the so called "p- value" of C/f • The "p- value" is defined 
as the probability that the relevant statistic takes a value 
at least as extreme as the one observed by the data when 
the null hypothesis is true. We have found for C7| that the 
"p-value" is 15.0%. 

Then we can conclude that even if we have a relatively 
strong Ul at the lowest signal-to- noise ratio, it is not un- 
likely to have such a high value by chance. Therefore, even 
considering the results from the coaddition map-making, 
Archeops data is still compatible with our Gaussian simu- 
lations. 

Although the high value found for C/f for the coaddi- 
tion map is not significant enough to be incompatible with 
Gaussianity, it is clear that there is a steady increase of U% 
when s/n decreases. This suggests the presence of system- 
atics in the coaddition maps that can depend on the resolu- 

3 rf(y)dy 
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Fig. 9. U 2 statistics of Mirage Archeops Data for different 
signal-to-noise cuts. 




Fig. 10. U 2 statistics of coaddition Archeops Data for dif- 
ferent signal-to-noise cuts. 

Table 3. U 2 from Archeops Mirage (coaddition) map for 
(s/n) c = 0.30 ((s/n) c = 0.27) and the probability that 
one s + n Gaussian Mirage (coaddition) simulation has a 
U 2 statistic larger than those of data. More precisely, the 
probability for J/f in the coaddition case is 0.6%. 
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tion. Moreover, the fact that it only appears in coaddition 
data suggests the possibility that it is a map-making issue. 
This also implies that systematics are better controlled in 
the Mirage than in the coaddition map-making. Therefore 
hereafter we focus only on the Mirage map making data. 

We performed a y 2 test with the three Minkowski func- 
tional using 11 thresholds from — 2.5<r to 2.5<r. We anal- 
ysed the Mirage data and a set of 1000 CMB Gaussian 
simulations with noise of the Mirage type. The correspond- 
ing histogram of the x 2 values of these simulations and of 
the data are presented in figure ITTi As it can be seen, the 
data are compatible with the Gaussian simulations. 

5.1. Systematic and foreground contamination 

The R&BT can also provide a powerful tool for estimating 
the level of this contribution. The test consists in adding 
different percentages of a template map to the Archeops 
143K03 bolometer map, for the Mirage and coaddition sim- 
ulations cases, to compare the resulting U 2 statistics to the 
ones obtained with the Archeops data at 143 GHz. 

This template map is c omputed from t he co added 
Archeops 353 GHz map (see iPonthieu et all |2005). This 
map contains thermal dust emission, atmospheric residuals 




10 20 30 40 50 60 70 



Fig. 11. Distribution of the x 2 values from the Minkonwski 
Gaussianity test for Archeops Mirage map. Vertical line 
shows the data results. Their cumulative probability is 
83.9%. 




Fig. 12. Mean of 10 4 Z7f statistics, from 10 4 signal plus 
noise Mirage simulations plus a factor times the con- 
tamination template. 0.0 < < 0.5. 

as the dominant components and also instrumental noise 
and CMB residuals. Thus, extrapolated to 143 GHz it will 
provide a good template of what could be a dust plus at- 
mospheric contamination at this frequency. 

Thermal dust is assumed to have a grey-body 
emission: v 2 B{y) which can be appro ximated in the 
Rayleigh- Jeans domain to Trj oc v 2 (see IPonthieu et all 
2005). Atmospheric residuals emission law has been es- 
timated empirically by the Archeops collaboration (see 
iMacfas-Perez et all 120071 ) and is also proportional to v 2 in 
the Rayleigh- Jeans domain. Dust and atmospheric resid- 
uals being the two main components, Archeops 353 GHz 
map has been extrapolated to 143 GHz by assuming that 
emission power law. Due to the extrapolation the CMB con- 
tribution on the 353 GHz template map is negligible with 
respect to the CMB at 143 GHz. 

U 2 statistic is the most sensitive to this effect as can 
be seen in figure [12] for the Mirage case where this statis- 
tic presents a prominent peak at signal-to-noise ratio cuts 
around 1.88. 

In order to determine the level of contamination we 
performed a x 2 test with the U2 statistic computed at 
(s/n) c = 1.88. It is optimal to perform a x 2 test with U2 
because U2 is normally distributed for the null hypothesis. 
Thus we can define 

X 2 (a d ) = - T ±—(U 2 -(U 2 )* d ) 2 (12) 

where (U2)a d and cr ad (U2) are the mean and the dispersion 
of U2 for CMB Gaussian simulations with noise plus a factor 
ad times the contamination template. In the left panel of 
figure [13] we present the x 2 of Archeops Mirage data for 
different 01.&. We can see that the minimum x 2 (best fit) 
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Fig. 13. From left to right, x 2 value of Archeops data for different ad and the histogram of best fit ad for a set of 
1000 Gaussian simulations without dust. These results have been obtained with the U2 statistic. 



occurs for a^=0.0. Analysing Gaussian simulations without 
dust we find that most of them reach the best fit for low 
values of ad (right panel of figure [T3|) . Specifically we have 
that a d < 0.27 for 90% confidence level (CL), and a d < 0.33 
for 95% CL. By comparing the dispersion of both maps, 
Archeops and 0.27 times the contamination template, we 
can exclude a dust plus atmospheric contamination larger 
than 7.8%. 

We computed another x 2 statistic using the Minkowski 
functionals for the dust analysis. In this case 

X 2 (a d ) = £(„(i) - {v(i)) ad )C^(vU) ~ (v(j)) ad ) (13) 



i and j cover 11 thresholds from —2.5a to 2.5a and the 
three Minkowski functionals. (v(i)) ad is the mean value of 
the corresponding functional at the corresponding thresh- 
old for Gaussian CMB simulations with noise plus ad times 
the dust template. C is the covariance matrix for Gaussian 
CMB simulations with noise. The value of ad that best fits 
Archeops data is ad = 0.0. Analysing Gaussian simulations 
without dust we find that ad < 0.28 for 90% CL, and ad < 
0.35 for 95% CL. 



5.2. Primordial non-Gaussianity 

There are several possible inflationary scenarios in 
which the primordial fluctuations are not Gaussian dis- 
tributed. The idea is to work with a simple non- 
Gaussianity model and to impose some constraints on 
it. In parti cular, we consider the "weak non-linear cou- 
pling; case" dKomatsu fc Sperget 1200 ll : iLiguori et all 120031 : 
iBartolo et all l2004f ) 



= $ L (x) + fni{$ 2 (x) - ($ 2 (x))} 



(14) 



where <fr(x) is the primordial gravitational potential, (which 
satisfies (4>(x)) = 0), $>l(&) is the linear random compo- 
nent (Gaussian distributed), and f n \ is the non-linear di- 
mensionlesfl coupling parameter. 

Scales larger than 1 degree are larger than the horizon 
scale at the recombin ation time, when CMB was formed 
([Liddle fc Lvthl |2000[ ). In this regime it is possible to 
make a good approximation linking CMB fluctuations and 
ravitational flu ctuations through the Sachs- Wolfe effect 
Sachs fc WolfcL I1967D AT(n)/T = $(n)/3 (notice how- 
ever that a better approximation should include the inte- 
grated Sachs- Wolfe effect). 



We analysed signal plus noise simulations with a f n \ 
term in this way, 

AT» = AT s (n) + ^{AT s (n) 2 - (AT s (n) 2 )} 



AT(n) = AT s '(n) + AT n (n) 



(15) 



where AT S is a Gaussian signal simulation, AT n is a 
Gaussian noise simulation, T = 2.725 K and AT is the 
analysed simulation. 

We performed a x 2 analysis for the primordial non- 
Gaussianity similar to the dust case for both /7 2 and the 
Minkowski functionals. The signal-to- noise eigenmodes yi 
are weakly dependent on f n \ . It can be seen that the mean 
value of yf for simulations with f n \ is 

(Vi)fni = 1 + 7—775 S * f nl + 1 1 m S * & ( 16 ) 

I + {JJA)i I + {VA)i 

*i = Y Y,^ R A L N 1 )ij^A) + (SksfyiLrfRA)* (17) 



(18) 



3,k 



4 We use the units system with c — 1. 



where bi is about an order of magnitude larger than cli for 
most of the s/n eigenmodes. This implies that (yf) f nl — l~ 
O(fni) which explains the low sensitivity of U2 to f n \ 
variations. In particular, we have found that it is much 
less sensitive than the Minkowski functionals. If we con- 
sider for example, a value of f n \ = 2300, we find a rel- 
ative variation ((yf)f nl — (y 2 )o)/(y 2 )o — 0.05 (and there- 
fore a similar ratio for U2 and E/f) for the former and 
(( f2 )Ui ~ ( f2 )o)/(F 2 ) ~ 0.50 for the latter. 

Therefore we performed a x 2 test with the three 
Minkowski functionals using different thresholds between 
— 2.5<j and 2.5<r. In the left panel of figure [TH we present 
the x 2 value of the data for different f n \ cases. We can 
see that the minimum x 2 value is reached for f n \ = 
200. Taking also into account the results obtained when 
analysing Gaussian simulations (see right panel of figure 
[14]) we can put the following constraints on f n \ from the 
Archeops data: f nl = 200 ± 6 3 ° ° at 68% CL, f nl = 200 ±f ° 
at 90% CL, and f nl = 200 ±|Jg° at 95% CL. 

6. Complementary analysis: WMAP in the same 
region 

WMAP is a NASA satellite dedicated to observe the 
anisotropics of the CMB with high accuracy at five differ- 
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Fig. 14. From left to right, x 2 value of Archeops data for different f n \ and the histogram of best fit f n \ for a set of 
1000 CMB Gaussian simulations with noise. These results have been obtained with the Minkowski functionals. 



ent frequencies between 23 and 94 GHz. Scientific results 
of this mission have allowed us to have a clearer image of 
the early universe, and to reduce the uncertainties in sev- 
eral cosmological parameters. Data products of this mission 
can be found on the wet0. 



6.1. The WMAP data 

We have analysed WMAP data with the same goodness- 
of-fit and the Minkowski functionals tests already used 
on Archeops data. The main purpose of this analysis is 
to compare Archeops results with a different experiment 
to discriminate among systematics, foreground emissions 
and intrinsic CMB no n- Gaussian features. It is clear that 
the WMAP frequencies complement very well the Archeops 
ones. A detailed analysis of the possible WMAP non- 
Gaussianities with this goodness-of-fit method deserves an- 
other work. 

The maps that we have analysed have been produced 
from the 1-year and 3-year WMAP foreground cleaned 
maps for the differencing assemblies corresponding to 
the cosmological frequencies 40, 60 and 90 GHz. The 
main properti e s of th e se maps are describe d in detail in 
iBennett et all (|2003aD : iHinshaw et~aTl (|2007f ) respectively. 

Specific ally we have u s ed the "combined map" as de- 
scrib ed in I Bennett et al.l ([2003a), (see also IVielva et al.L 
I2004D . The WMAP CMB simulations which are used in 
the analysis are also combined simulations, that is, CMB 
signal simulations were produced for each channel and then 
combined in the same manner than for the data. 

According to IBennett et~aD (|2003aD WMAP noise is 
highly uncorrelated, that is, the noise from a given pixel 
i is independent of the noise from another pixel j. The 
noise combined simulations are produced fro m the "com- 
bined variance map" as it is shown in e.g. IVielva et al.l 

dlooi). 

We have analysed both combined maps, 1-year and 3- 
year (hereafter WCM1 and WCM3). The WMAP mask 
considered for both analyses was the 3-year KpO one be- 
cause it is the most conservative on e for WCM3 a n d also 
contains the 1-year KpO mask. See Hinsha w et al.l (|2007f) 
for details about new masks and IBennett et al.l (|2QQ3bh for 
original masks. The actual mask we have used is the 3-year 
WMAP KpO degraded to our resolution times the Archeops 
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Fig. 15. WCM3 Data at HEALPix resolution N side = 32 
(it corresponds to a pixel size of ~ 1.8 degrees). This map is 
centered on Galactic longitude / = 180 degrees. The pixels 
contaminated by Galactic and extragalactic emission are 
covered with the mask described in the text. Grid lines are 
spaced by 20 degrees. 



mask0. Its number of pixels is 1648. In figure [15] WCM 3 data 
is plotted using this mask. 

6.2. Gaussianity test on WMAP data 

In order to perform the R&BT test on WCM1 and WCM3 
maps we followed the same steps than for the Archeops 
analysis. We calculated their corresponding S and N matri- 
ces for the 1648 pixels available after applying the combined 
Archeops- WMAP mask. 

We assume the best fit model of the 3-year WMAP data 
for both analysis, WCM1 and WCM3. At the resolution 
with which we are dealing, 1.8 degrees, the power spectra 
of the 1-year and 3-year data are very approximately the 
same. This assumption implies that the S matrix is the 
same for both releases. The S matrix is computed from 
1.2 x 10 5 Gaussian simulations following equation [6l Each 
simulation was produced in the same 90 dual core proces- 
sors mentioned before, and took an average CPU time of 
360 s and an average RAM memory of 0.4 GB. 

As commented above, WMAP noise is highly uncorre- 
lated and therefore we can assume that the noise matrices 
are diagonal. This means that the correlation element cor- 
responding to pixels i and j is Nij = af * Sij, where af 
is the combined noise of pixel i. Noise matrices for WCM1 



http : //lambda . gsf c . nasa . gov/ 



6 For comparison, we have also repeated the goodness-of-fit 
analysis on Archeops data using this combined mask, finding 
similar results to those obtained in section 5 using the Archeops 
mask. 
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Fig. 16. From left to right, mean and dispersion of Uf statistics (where i goes from 1 to 5) for different signal-to-noise 
cuts, corresponding to 10 3 signal plus noise WCM3 simulations. 
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Fig. 17. From left to right, Uf statistics for WCM1 and WCM3 presented for different signal-to-noise cuts. 



and WCM3 must be constructed with their corresponding 
noise variances which differ by an approximate factor of 3. 

Two additional sets of 10 3 Gaussian signal plus noise 
simulations (corresponding to WCM1 and WCM3 maps) 
were performed for the calibration of the matrices. In fig- 
ure HH we present the mean and the dispersion of the Uf 
statistics at different signal-to-noise cuts for the WCM3 
case. Note that the numerical range for the possible signal- 
to-noise cuts (s/n) c is wider than for the Archeops case, 
because WCM3 noise is smaller than the Archeops one at 
this resolution. (s/n) c range for WCM1 is approximately 
the same than that of WCM3 reduced by a factor y/3. The 
mean and the dispersion for WCM1 simulations are simi- 
lar to those obtained for WCM3. It can be seen that mean 



Table 4. Mean and dispersion of Uf statistics from 10 3 
WCM1 simulations for (s/ri) c = 3.64. 
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Table 5. Mean and dispersion of Uf statistics from 10 3 
WCM3 simulations for (s/n) c = 6.33. 
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values of Uf statistics are close to one for almost all signal- 
to-noise cuts and all the computed statistics, but the dis- 
persion becomes higher than square root of two for high 
signal-to-noise cuts and for statistics with high order mo- 
ments, like E/f and higher order statistics. 




Fig. 18. Distribution of the x 2 values from the Minkonwski 
Gaussian test for WCM3 data. Vertical line shows the data 
results. Their cumulative probability is 12.0%. 



As for the Archeops case, these high values are explained 
by the small errors present in the computed correlation ma- 
trices plus small numerical errors in the diagonalization of 
these matrices, which are amplified through the high order 
moments. In table 2] we present the mean and the disper- 
sion of Uf statistics for 10 3 WCM1 simulations with noise 
for all the eigenmodes (s/n > 3.64). Note how the disper- 
sion is increasing with the order of the statistics. In table [5] 
the same quantities are presented for 10 3 WCM3, obtained 
also from all the eigenmodes (s/n > 6.33). The effect is 
the same for the high order moment statistics. The results 
for the Uf statistics for WCM1 and WCM3 data maps are 
presented in figure [171 As can be seen, all Uf values satisfy 
Uf < 7.15. The upper limit 7.15 corresponds to a upper 
tail probability of 0.7% for the theoretical distribution. In 
order to confirm or rule out a possible non-Gaussian de- 
tection, this result should be studied more carefully. First 
of all, we have that for both WCM1 and WCM3 U% is the 
only statistic which reaches some sharp peaks above 6.6 
(which corresponds to a upper tail probability for the the- 
oretical distribution of 1.0%). From the plots in figure [T7J 
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Table 6. WCM1 Uf statistics for (s/n) c = 21.81, and their 
corresponding upper tail probabilities. 
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Table 7. WCM3 Uf statistics for (s/n) c = 37.92, and their 
corresponding upper tail probabilities. 
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Ul reaches this peak at (s/n) c = 21.81 for WCM1 and 
(s/n) c = 37.92 for WCM3. We estimated the upper tail 
probability for the Uf statistics of data at the mentioned 
signal-to-noise cut by performing 10 3 Gaussian simulations. 
These results are presented in tables [6] and As we can see 
for the Ul statistic, we have that this probability is 1.0% 
and 0.7% for WCM1 and WCM3 respectively, very similar 
to the theoretical value. 

This probability is obtained for the precise signal-to- 
noise cut where f7f reaches its maximum. Since the width 
of the maxima is much smaller than the range of variation 
of the signal-to-noise eigenvalues, it makes sense to ask for 
the significance of the detection. Thus, from the simulations 
we computed the "p- value", i.e. the probability of finding 
a value of f7f larger than 7.15 at any signal-to-noise cut, 
the maximum value reached by the data. This probability 
is 18% for WCM1 and 17% for WCM3. 

From the previous discussion, we conclude that the 
sharp peaks found in the data are not significant. Also, well 
studied cases of artificial CMB non-Gaussianities, like skew- 
ness or kurtosis produced using; the Ed geworth expansion 
(see iMartmez-Gonzalez E. et "all (|2002f ) for applications of 
this expansion to the CMB no n- Gaussianity analyses), usu- 
ally show deviations of the Uf statistics in the form of a 
large plateau. Besides, we would like to remark that at the 
signal-to-noise cuts where the maxima are found there are 
less than one hundred {yi} numbers to compute the Uf 
statistics (around 70), and the test works correctly only 
asymptotically (n » 1). 

WCM3 data were also analysed with the Minkowski 
functionals as in the Archeops case (that is, using 11 thresh- 
olds between — 2.5<r and 2.5<r and the three functionals). 
The histogram corresponding to the x 2 values for 1000 
Gaussian simulations and the value for WCM3 data are 
presented in the figure [TBI As we can see the WCM3 data 
are again compatible with Gaussianity. 

Finally, we performed an analysis on simulations with 
f n i parameter as defined in equation [151 The procedure 
was the same as the one performed for Archeops case. As 
discussed in section 5.2, we only use the Minkowski func- 
tionals for the f n i case. The x 2 value for WCM3 data is 
minimum for f n \ = 100. Analysing Gaussian simulations, 
the constraints found for j n \ are: f n i = 100 ^200 a ^ 68% 
CL, f nl = 100 tfo at 90% CL, f nl = 100 at 95% 
CL. These limits are compatible with those obtained from 
Archeops since the tighter constraints found for WCM3 can 
be explained by the significantly smaller noise in that ex- 
periment. In particular, if we analyse simulated Archeops 



data with noise normalized to the same amplitude as that 
of WCM3 we find similar limits for j n \. 

7. Conclusions 

The expected behaviour of the Uf statistics as a Xi dis- 
tribution has been confirmed for the order index interval 
1 < i < 4 with "realistic" simulations assuming Gaussian 
CMB anisotropics. For higher moments, i > 4, the mean of 
the distribution is \i ~ 1 but the variance is a > 2. This 
is because of the propagation of errors through higher or- 
der moments which in practice complicates the use of high 
order Uf in our analysis. 

From the analysis of both kind of Archeops maps, coad- 
dition and Mirage, we have found that both are compatible 
with Gaussianity. Only the f7f statistic for coaddition map 
is close to 8.0 for low (s/n) c . Although in principle the 
probability that Ul takes values bigger than 8.0 for a given 
signal-to-noise cut in the Gaussian hypothesis is very low 
(see table [3]), the corresponding "p- value" for having Ul 
larger than 8.0 at any signal-to- noise cut is 0.1482. This 
is not negligible and thus this detection is not significant. 
Moreover this effect does not appear in the Mirage map, 
and therefore should be assigned to issues related to the 
map-making process. 

The analysis with the Minkowski functionals on the 
Mirage map also returns compatibility with Gaussianity. 

Our analysis also implies constraints on the amount of 
contamination that can be present at 143 GHz. Using as 
template for dust and atmosphere the Archeops map at 
353 GHz, we limit the possible contamination to be lower 
than 7.8% at 90% CL using U2 statistic. A similar limit 
is obtained with the Minkowski functionals. 

We have also compared the Archeops results with the 
WMAP 1 and 3-year data in the same region of the sky. 
For both sets of data a sharp peak in U% has been found 
at specific signal-to- noise cuts. Although the probability of 
finding such a peak at a given signal-to-noise cut is very 
small, the "p-value" obtained when different cuts are al- 
lowed is appreciable. Therefore we can conclude that the 
WMAP data, when the same region than Archeops is con- 
sidered, are also consistent with Gaussianity. The same 
conclusion is reached when the data are analysed with the 
Minkowski functionals. 

Finally, we have established a constraint in the value of 
the non-linear coupling parameter f n \ . Analysing Archeops 
data, we found that f nl = 200 ±f ° at 90% CL, and f ni = 
200 

-800° a ^ 95% CL. When the same analysis was done 
with WMC3 data using Archeops- WMAP combined mask, 
we found f n i = 100 ±f ° at 90% CL, f nl = 100 at 95% 
CL. These limits are similar to the ones expected for an 
Archeops-like experiment with a noise amplitude similar to 
that of WCM3. 
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We have used the CAMB code ([Lewis et all 120001) for our analysis . 
The CAMB code is derived from CMBFAST (|Zaldarriaga & S eliak, 
2000 ). The HEALP ix package was used throughout the data analysis 
(|G6rski et all 120051) . 
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